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ABSTRACT 

We study relativistic particles undergoing surfing acceleration at perpendic- 
ular shocks. We assume that particles undergo diffusion in the component of 
momentum perpendicular to the shock plane due to moderate fluctuations in the 
shock electric and magnetic fields. We show that dN/dE, the number of surfing- 
accelerated particles per unit energy, attains a power-law form, dN/dE oc E~^. 
We calculate b analytically in the limit of weak momentum diffusion, and use 
Monte Carlo test-particle calculations to evaluate b in the weak, moderate, and 
strong momentum-diffusion limits. 

1. Introduction 

The acceleration of high-energy particles is an important problem in astrophysics. This 
paper examines one acceleration mechanism, "surfing acceleration,"^ as it applies to relativis- 
tic particles at perpendicular shocks. An illustration of an ion undergoing surfing acceleration 
at a perpendicular shock is given in figure 1. In this picture an ion of velocity v arrives just 
upstream of the shock with small \vx\ and is reflected by the jump in the electrostatic po- 
tential $(a;) at the shock, which is illustrated schematically in figure 2. The Lorentz force 
then causes the ion to return to the shock where it is again reflected and brought back to 
the shock front by the Lorentz force. This process continues, confining the particle to the 
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et al (2001), Lever, Quest, & Shapiro (2001), Shapiro, Lee, & Quest (2001), and Hoshino & Shimada (2002). 
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vicinity of the shock, where it is accelerated in the — y direction by the motional electric 
field, — u X B/c, where u is the plasma velocity. If one thinks of q^{x) in figure 2 as analo- 
gous to the profile of an ocean wave, the accelerating ion in figure 1 is analogous to a surfer. 
In this paper, we will focus on ions. However, shock surfing acceleration is believed to be 
important for electrons as well (McClements et al 2001, Hoshino & Shimada 2002), and it 
is trivial to modify the analysis of sections 2 and 3 to treat the electron case. 
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Fig. 1. — Idealized trajectory of a surfing ion at a perpendicular shock as seen in the shock 
wave frame. 

An ion's x-oscillations in front of the shock can be approximately described in terms of 
an effective potential energy (Lee et al 1996, Ucer & Shapiro 2001), which can be derived as 
follows . We start with the x-component of the equation of motion, 

„| (,.,^^M.). ,) 

where 7 is the relativistic Lorentz factor, and wc work in shock wave frame. We assume 
that Ex = —d^/dx with $ a function of x alone, that l^z-Bj/l <^ \vyBz\^ that \vx\ <^ and 
that the particle oscillates in x on a time scale that is much shorter than the time required 
for 7 or Vy to change appreciably. We also treat as spatially uniform. We then multiply 
equation (1) by Vx and integrate in time, treating 7 and Vy as constant, to obtain 

2 

Tfi'yv 

— — + Ues{x) — constant =w-\- C/mm, (2) 
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Fig. 2. — Idealized plot of electrostatic potential energy q^{x) and effective potential en- 
ergy U^six). 



where 



Ues{x) = q 



'. , XVyB^ 

— 



(3) 



w is the energy of x-oscillations, and C/min is the value of Ues at its local minimum at x — x^ 
We define 



(4) 



where f/max is the value of Ues{x) at its local maximum. Surfing ions, which are trapped in 
the effective potential well, have < w < Wmax- The effective potential energy is plotted in 
figure 2 for Vy < 0. The turning points of an ion of given w, obtained by solving equation (2) 
with Vx — 0, are denoted xl{w) and xji{w). Over many bounces, 7 and Vy vary slowly, 
causing w to change in a way that can be calculated by noting that the adiabatic invariant 

/•xr{w) 

J = / p^dx (5) 

is approximately conserved (Ucer & Shapiro 2001). 

Surfing can only occur if there is a local minimum in Ues, which requires that 



VyB^ 



(6) 
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over some interval of x [see, e.g., Sagdeev & Shapiro (1973), Katsouleas & Dawson (1983), 
Ucer & Shapiro (2001)]. In order for particles with \vy\ ~ c to surf, \Ex\ must locally 
exceed B^, a condition that may arise at large Alfvenic Mach number (Hoshino & Shimada 
2002). 

The escape of a surfing particle from the effective potential well determines how far the 
particle can travel along the — y direction, which in turn determines the amount of energy 
the particle gains. Several escape mechanisms have been considered in the literature. First, if 
the maximum value of E^, denoted E^^ max) is less than (assumed constant) and a surfing 
particle's \vy\ increases to the point that \vy\ > cEx^max/Bz, then the local minimum of Ucs 
disappears and the particle escapes downstream (Sagdeev & Shapiro 1973). Second, as a 
non-relativistic particle's \vy\ increases, the slope of Ues at large x increases as illustrated in 
figure 2, causing the width Ax of a particle's bounce oscillations to decrease. Since J ~ Axp^ 
is approximately conserved, where is a typical value of \px\, decreasing Ax increases p^, 
thereby increasing w (Ucer & Shapiro 2001). If w increases above Wmax, then the particle 
escapes downstream (Lee et al 1996, Ucer & Shapiro 2001). Third, in an oblique shock the 
velocity vector of a surfing particle gradually rotates in the yz-plane due to the Lorentz 
force associated with B^, as depicted in figure 3, which is an adaptation of figure 4 of Lee 
et al (1996). As this happens, \vy\ is eventually reduced to zero and the particle escapes 
upstream since the VyB-^/c component of the Lorentz force ceases to turn the particle back 
towards the shock (Lee et al 1996). Fourth, wave-particle interactions can reduce a surfing 
particle's \vy\ to zero, again allowing the particle to escape upstream (Shapiro, Lee, & Quest 
2001). Fifth, as noted by Lee, Shapiro, & Sagdeev (1996) and Hoshino & Shimada (2002), 
some shocks are intrinsically nonstationary and can periodically "break," allowing surfing 
particles to escape downstream. 

For perpendicular shocks and relativistic particles with \vy\ ~ c, the first three mech- 
anisms described above are absent. This has led previous authors to consider scenarios of 
"unlimited acceleration" when E^ > B^ over some interval in x (Katsouleas & Dawson 1983, 
Ucer & Shapiro 2001). 

The purpose of this paper is to show that shock surfing of relativistic particles at per- 
pendicular shocks leads to a power-law energy spectrum of accelerated particles, provided 
fiuctuations in the shock electric and magnetic fields are not too large, and to calculate the 
power-law index of this spectrum. We ignore diffusion in a particle's ^/-momentum, py, and 
assume that changes in Py are dominated by the motional electric field. However, we take 
into account moderate fiuctuations in the electric and magnetic fields in the shock's vicinity, 
and in particular the effects of these fluctuations on the x component of a particle's momen- 
tum, px- An ion propagating primarily in the — y direction encounters a series of random 
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forces in the x direction (denoted 5 Fx) due to fluctuations in and B^i as illustrated in 
figure 4. These forces induce a series of random increments in p^., causing diffusion in p^. 
For example, suppose the typical electric and/or magnetic field structures have a length ly 

along the y direction and a typical amplitude ^Ej^, that an ion moves primarily in the — y 
direction at a speed ~ c, and that these structures fluctuate on a time scale that is longer 
than the time 

At ~ ^ (7) 
c 

for a surfing particle to traverse one of the structures. The random increment in when an 
ion moves through one structure then has a magnitude 

Ap^ ~ AtqSE^. (8) 

If the field structures are space filling, the ion undergoes one such momentum kick during 
each time At, leading to a coefficient of diffusion in p-^, 

Dp^ ^ > (9) 

that is independent of particle energy. The energy spectrum of accelerated particles can be 
understood in terms of two competing effects. On the one hand, diffusion in p^ stochas- 
tically de-traps particles by increasing w above Wmax, thereby allowing particles to escape 
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Fig. 4. — Simple visualization of the origin of diffusion in p^. 



downstream. On the other hand, the gradual increase in 7 due to acceleration by Ey focuses 
a surfing particle with l^^^l ~ c into the bottom of the effective potential well, as follows from 
conservation of J (Ucer & Shapiro 2001). This combination of factors leads to a power-law 
energy spectrum of accelerated particles, dN/d^ oc 7"''. The value of b depends upon 



C ^acc 

a — ~ , 



(10) 



where tacc = i^c'-f /q\Ey \ is the time to double a surfing particle's energy and tjig = m'-fw^s^x/ Dp^ 
is the momentum diffusion time scale. As a is decreased, the acceleration time is reduced 
relative to the momentum diffusion time scale, particles experience larger energy increases be- 
fore escaping the surfing mechanism, and one expects the energy spectrum dN/ 0^7 to become 
harder. This expectation is born out by the calculations and simulations to be presented 
in sections 2 and 3, which show that in the strong-scattering limit (a ^ 1) 6 ~ a, in the 
moderate-scattering hmit (a ~ 1) 6 ~ 1 — 3, and in the weak-scattering hmit (a <^ 1) 6 ~ 1. 

If, contrary to what we assume, the field fluctuations were so large that the minimum 
of the effective potential energy disappeared at locations separated by a typical distance Ay 
along a surfing particle's trajectory, one would not obtain a simple power-law energy spec- 
trum; the spectrum would instead steepen above an energy ~ q\Ey\Ay corresponding to the 
energy a particle gains by surfing a distance Ay (Hoshino & Shimada 2002). 
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As mentioned previously, if a particle is assumed to undergo diffusion in py with diffusion 
coefficient Dp^ as well as "advection" in py due to Ey, then a particle can in principle diffuse 
to Py = Q and escape into the upstream region. The condition under which this escape 
mechanism can be neglected can be understood by considering an analogy to first-order 
Fermi acceleration at a non-relativistic shock, as follows. Suppose a particle is a distance x 
downstream of a shock, that the downstream fluid velocity is Ux, and that the particle diffuses 
in space with diffusion coefficient D^. The time for a particle to return to the shock in the 
absence of advection is tdifr,x ~ / and the time for the particle's distance from the shock 
to be doubled by advection in the absence of diffusion is iadv — x/ux- The probability that the 
particle returns to the shock is e~*diff,x/tadv — ^-xu^/d^ ^ggg^ e.g., Bicout 1997). This advection- 
diffusion problem is analogous to the advection and diffusion in py that would occur in shock 
surfing with nonzero Dp^. The probability that a particle escapes the surfing mechanism by 
diffusing in momentum space to Py = Q is again ~ e-*diff,py/tadv,py ^ where tdifr.p,, = P^/^Py is 
the time to diffuse to Py = in the absence of advection in py, and tadv.p^ = ^acc = nT'C'y/q\Ey\ 
is the time for a particle's |/-momentum to double due to Ey. Escape through py diffusion is 
thus negligible provided iacc < tdis,py, or Dp^ < pf^qlEyl/mcy ~ q\PyEy\. 

We note that it is not clear how efficiently particles are injected into the surfing ac- 
celeration mechanism. Injection appears to be natural for pickup ions accelerated by the 
heliospheric termination shock. In the frame of the solar wind, pickup ions are born with 
particle speeds comparable to the solar- wind velocity, and thus pickup ions that arrive at the 
heliospheric termination shock at the right point in their gyromotion have very small \vx\, 
which allows them to be trapped in front of the shock (Lee et al 1996, Zank et al 1996, Lever 
et al 2001). Surfing acceleration in other cases, however, may be more difficult to initiate. 
It may be that particles can be injected at the time that a shock is formed if the shock 
originates from an explosive event. It is also possible that turbulent conditions in the shock 
transition region can feed particles into the surfing mechanism. Hoshino & Shimada (2002) 
have found electron surfing in PIC simulations of perpendicular shocks, demonstrating that 
injection does occur, at least for electrons. Additional work on injection, however, is needed. 

Although our analysis is restricted to the case of perpendicular shocks, some of the 
results may apply to the more general case of superluminal shocks, for which the point of 
intersection between a field fine and the shock plane moves faster than c. For such shocks, 
there exists a "perpendicular shock frame" in which the shock is stationary and the upstream 
electric and magnetic fields he in the shock plane (Begelman & Kirk 1990). There are, 
however, two issues that would have to be addressed before our results could be generalized 
to this case. First, the condition for E^ to exceed in more general superluminal shocks 
must be found. Second, in general the upstream fiow velocity in the perpendicular shock 
frame has a component along the upstream magnetic field. If particles start to surf with 
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a significant z velocity, they would only approach the state that we analyze with v oc — y 
after significant acceleration in the — y direction. As a particle approaches this state, \vy\ 
gradually increases, and the effects of such a gradual increase are not included in our analysis. 

The remainder of this paper is organized as follows. In section 2 we show analytically 
that surfing ions that diffuse in p^ are described by a power-law energy spectrum. We 
calculate analytically the power-law index for the case of weak momentum diffusion (a ^ 1), 
and derive a rough estimate for the power-law index in the case of strong momentum diffusion 
(a ^ 1). In section 3 we present Monte Carlo numerical calculations of surfing in the 
presence of diffusion in "p^ for the weak, moderate, and strong momentum-diffusion cases. 
We summarize our results in section 4. 



2. Analytic calculation of energy spectrum of accelerated ions 

We begin with the relativistic Vlasov equation with a term added to model diffusion 
in Vx, 



where D^^ is taken to be independent of particle energy as in equation (9). We consider a 
perpendicular shock and ignore spatial variations in the magnetic field and motional electric 
field: 

B = S,z, (12) 
Bz — constant > 0, (13) 



and 



We also take 



and 



Ey = constant < 0. (14) 



= 0, (15) 
c-^^^^^ = 0. (17) 



dt dy dz 

We confine our analysis to relativistic surfing ions, for which we assume 



Vy ~ -C, (18) 

\vx\ < c, (19) 
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Vz = 0, 
dy dz 

VxBz 



and 



7 = 



0, 

\Py\ 
mc 



Equation (11) can then be rewritten as 



Px 



dpi' 



We now change variables, defining 



t' = t, 



X = X, 



and 



or, equivalently. 



where 



m'yv'i TT f \ 



W — 



+ t/eff(a;) - t/min, 



(20) 
(21) 

(22) 
(23) 

(24) 



(25) 

(26) 
(27) 

(28) 
(29) 
(30) 



is the effective potential energy for the particle's motion in x illustrated in figure 2, w is 
the energy associated with this motion, and C/min is the value of U^s at its local minimum 
at x = Xmin- As discussed in the introduction, surfing particles satisfy Q < w < tfmax, where 



(31) 



and C/max is the value of U^s at its local maximum. In terms of (i', x\ w,Pyo), equation (24) 
becomes 



9/± , q\Ey\vldf^ d_( df 



± 



dt' dx' 2 c dw dw \ dw 

where /"*" (/~) is the distribution function for particles with > (< 0). 



(32) 
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We assume that the time tb for a particle to bounce in the effective potential well is 
much shorter than the acceleration time or the time for a particle to escape from the potential 
well due to diffusion in p^. We thus take 



and 



with 



We then set 



e= — <1. (35) 
To 

f^-fo+ef^ + --- (36) 

and expand equation (32) in powers of e. Upon collecting all terms of order e^/o/rg (there 
is only one) and dividing by \vx\, we find that 



dx' 



= 0. (37) 



As in the introduction, we take xl{w) and xr{w) to be the bounce points of a particle 
of "energy" w, with x^ < xj^. Since f^{xL{w),w) and f''{xL{w),w) give the distribution 
function at the same point in phase space, one has /+(a;/,(u'), w)) = {xl{w),w). Similarly, 
f^{xR{w),w) = f~'{xR{w),w). Equation (37) thus implies that 

/o+ = U = fo. (38) 

Upon collecting all terms in equation (32) of order e/o/rB and dividing by \vx\, we find that 

I dfo dfl _ q\Ey\\v^\ dfo _ ^ d_f, M\ 



\vx\ dt' dx' 2c dw dw \ dw ^ 

We integrate equation (39) from xl{w) to xr{w) for /+, then integrate equation (39) from 
xl{w) to xr{w) for /~, and then add the two resulting equations to annihilate the terms 
involving thereby obtaining 



9/o q\Ey\g dfo _ ^ d f [2 dfo 



dt' c^j2rrvy dw dw \\m'y dw ' ' ^^^"^ 

where 



9=1 dx \Jw- Uesix) + C/min, (41) 
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and 



tb 



I 



dg 
dw 



(42) 



'xl{w) Fa; I 

We note the behavior of g{w) near w = for future reference. Near x — x„im, 

Ues{x) = [/min + a{x - Xmin)^ + P{x - X^iinf + k{x - Xmin)^ + ■ ■ ■ (43) 

Using the technique of asymptotic matching, one can show that as w — > 0, g{w) has the 
asymptotic expansion 



g{w) = 



TTW 



(44) 



Since dg/dw exists and is integrable, the asymptotic expansion of dg/dw near w = is given 
by term-wise differentiation of equation (44) : 



TT 



dg ^ 

dw 2y/a 



(46) 



We now recast equation (40) into the standard Fokker-Planck form for the number of 
particles per unit w per unit Py^ , 



( f+ 



{f+ + f-)dx 



^ 2tb/o. 



xl{w) 



(46) 



We start by noting from equations (23), (25) and (27) that 

-Pyo - qEyt' 



7 



mc 



where both pyo and Ey are negative. Equations (42) and (47) give 

dTB _ TBq\Ey\ 



Next, we define 



and note that 



dt' 2'ymc 
w dg 



h 



g dw^ 



d f P 



dw \2tb J 2v^2m7 



d ( Pw 



dw V h 



-P 



(47) 

(48) 

(49) 
(50) 



- 12 - 



Throughout the rest of this section we drop the primes, setting t' ^ t. Substituting /o = 
P/2tb into equation (40) and making use of equations (48) and (50), we find 

^^"1^' ^^^^ 

where 

^_\^] p_l_ f Dv.^P \ (52) 
\ m7 2mc'-fh J dw \ m'-fh J 

is the flux of particles in w-space. The quantity in the first set of parentheses on the right- 
hand side of equation (52) is (Aw) /At, while the coefficient of P within the second set of 
parentheses is ((Aw)^)/2At, where Aw is the change in w during a small time interval At, 
and the angle brackets denote an average over the stochastic scattering that gives rise to dif- 
fusion inp-r. The term —q\Ey\w/2mc'jh is the rate of change of w caused by the approximate 
conservation of the adiabatic invariant 



J ^ Pxdx ^ g-s/2mrf. (53) 



That is, setting dJ/dt = yields dw/dt = —q\Ey\w/2mc^fh. which confirms the statement 
made in the introduction that J-conservation drives a relativistic surfing particle towards the 
bottom of the effective potential well (Ucer & Shapiro 2001). The terms on the right-hand 
side of equation (52) proportional to Dp^ give the advection and diffusion in w-space arising 
from a time average of a particle's diffusion in px- 

If diffusion in px causes a particle's w to increase to Wmax, the particle escapes down- 
stream. We thus seek a solution to equation (51) subject to the boundary condition 

P(«;max, t) = 0. (54) 

We also require that P be finite and differentiable aiw = 0. Since there are no PyQ derivatives 
in equation (51), we can treat particles with different values of PyQ independently. We will 
henceforth focus on particles with a single value of Pyo, and assume that a large number of 
such particles are injected into the surfing mechanism at t — 0. We rewrite equation (47), 
dropping the prime on t', to obtain 

7-f^)(rf + i), (55) 
\ mc J 

where d = pyo/qEy is the time required to double the particle's initial ^/-momentum. Be- 
cause a surfing particle's energy increases linearly in time, the number of particles that are 

dN 

accelerated to a final Lorentz factor in the interval (7,7 -|- ^7), denoted -j—dj, is equal to 

07 
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the number of particles that escape the effective potential well in the time interval {t,t + dt), 
where t — {mcy / q\Ey\) — d and dt — {■mc/q\Ey\)d'y from equation (55). Thus, 



dN ^ mcr(wmax, t) 
dl (l\Ey\ 



t={mc'y/q\Ey\)-d 



(56) 



Using equation (55), we can rewrite equation (51) as 



dP d 



dt du 

where 



u a 



&^ ( auP\ _ 



u = (58) 

^max 

and 

The eigenvalue equation 

C{Pn) = -XnPn (60) 

can be written in Sturm-Liouville form. We consider eigenfunctions that are finite and 
differentiable at m = and satisfy 

Pn{l) = 0. (61) 

For these boundary conditions, and since g{0) = 0, the eigenvalues are real and the eigen- 
functions of C form a complete set on the interval < m < 1. The solution for P{u,t) can 
thus be expanded as 

P{u,t) = J2Xn{t)Pn{u). (62) 

n=0 

Prom equation (57), 



Xn{t)=Xn{0)[j^^j , (63) 

where Xn(0) is determined from the initial conditions on P. We order the eigenvalues and 
their corresponding eigenfunctions so that Aq < Ai < A2 < . . . , etc. For t ^ d, P{u, t) is 
dominated by the first term in the sum in equation (62), 

P{u,t) (X PQ{u)t-^\ (64) 

and, from equations (52), (56), and (64), 

dN f, 

^cc,-\ (66) 
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where 



6 = 1 + Ao. 



We note that 



C{Pn) 



a d 
2du 



du \ g' 



where 



J _dg_ 
du 



9 = 



(66) 

(67) 
(68) 



Multiplying equation (60) by e'^^"Pn/g', integrating from u — to u — 1, and integrating by 
parts, we find that 

A„ > 0. (69) 



2.1. Smallest eigenvalue in the limit of weeik momentum diffusion (a <^ 1) 

We now calculate Aq to leading order in a when a <^ 1 for any Ues{x) with a single 
potential well as in figure 2. We take 

Po(0) = 1. (70) 

Prom Sturm- Liouville theory, we know that Aq is the unique value of A„ in equation (60) 
for which P„ satisfies the boundary conditions and has no zeroes in the interval < u < 1. 
We use this fact to calculate Aq in the following four steps: (1) we assume an order of 
magnitude for Aq, (2) we derive the resulting solution for Pq, (3) we determine the value 
of Aq by requiring that Po(0) = 1 and Po(l) = 0, and finally (4) we verify that Pq has no 
zeroes for < n < 1. 

Starting with step 1, we set 

Ao = a-^e"^/"^ (71) 

and assume that 6 is of order unity. This value is close to A's lower limit of 0. The correspond- 
ing particle energy spectrum is hard {b ~ 1), which is expected in the small-a limit, since 
the acceleration time is much less than the momentum-diffusion time [see equation (10)]. 

We now carry out step 2. We set 

^ (0) 

and rewrite >C(Po) = — AqPo as 

d_ 
du 



e""^"^^ (e"/"P) 



W^^-Va + g-l/a^)_ (73) 
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Equations (70) and (61) give 

F(0) = 0, (74) 

and 

F(l) = -1. (75) 

We assume that F decreases monotonically from to —1 os u increases from to 1, an 
assumption that we verify at the end. Then for |1 — ii| » a, we can take e~"/" + e~^/°F ~ 
Q-u/a -j-j^g right-hand side of equation (73). In addition, when |1 — w| <S 1, the entire 
right-hand side of equation (73) is of order e~^/" and can be replaced by to an excellent 
approximation. Thus, throughout the interval Q < u < 1, we can neglect the second term in 
parentheses on the right-hand side of equation (73), and write 



d_ 
du 



ou ^ ' 



(76) 



Integrating twice, using g{^) = from equation (44), and imposing equation (74), we find 
that 



/ du, —— / du2 e--'/-g'{u2). (77) 
Jo 9{ui) Jo 



We now carry out step 3. Imposing equation (75) and evaluating the integral in equa- 
tion (77) to lowest order in a, we find 

For step 4, we note that for < w < 1, and g' are positive. Thus, F decreases 
monotonically from to —1 as u is increased from to 1, as assumed. This implies that P 
has no zeroes for < rt < 1, and that the value of Aq given by equations (71) and (78) is 
indeed (a good approximation of) the smallest eigenvalue. 

Equations (45), (71), and (78) give 

where a is defined by equation (43) . The power-law index of the energy spectrum of accel- 
erated particles is then 

6=1 + y\_rv__ (^gQ^ 

TTWrnax 

Prom equation (72), it can be seen that Pq is strongly peaked near u — 0. This means that 
if a large number of particles start surfing at t — 0, the majority of those remaining at large 
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times will be concentrated near the bottom of the effective potential well (with u < a since 
Po ^ 6-"/°^ for u < 1). 

For reference, the higher eigenvalues can be found by expanding equation (60) in powers 
of a, solving for P/h in an inner region with u <^ 1 and an outer region with u a, and 
then matching the inner and outer solutions in the region a <^ u <^ 1. One finds that 
P/h satisfies Rummer's equation in the inner region regardless of the precise form of f/off (x). 
This is because /i ~ 1 in the inner region regardless of the exact form of Uef[{x). The inner 
solution can be matched to the outer solution only if 

An ^ I (81) 

for n = 1,2,..., provided n is of order unity (e.g.,^ o;"^). The n zeroes of P„, de- 
noted ui,U2, . . . ,Un, occur in the inner region for n of order unity, with each of the Ui 
of order a. 



2.2. Smallest eigenvalue in the limit of strong momentum diffusion (a ^ 1) 

When q; ^ 1, we can write equation (60) to lowest order in as 

^2 fuP\ d fP\ 

Stt - IT- o = -V'^, (82) 



du^ \ 2h J du\2 
where 

A = aij. (83) 

We again seek a solution for P{u) that satisfies P{0) — 1 and P(l) = 0. In going from 
equation (60) to equation (82) we have retained the highest order derivative. Wc thus have 
a regular perturbation problem and do not expect a boundary layer. Since there is no small 
parameter in equation (82), we expect the smallest value of ip to be of order unity. Then, to 
lowest order in a"^, 

h^a (84) 

to within a factor of order unity. 



3. Numerical calculation of energy spectrum of accelerated pairticles 



In this section, we confirm and extend the analytic results of section 2 with the use of 
Monte Carlo test-particle calculations. We follow a set of test particles (ions) with = 



that obey the equations of motion 

dpx 
dt 




+ 6p, 



'XI 



(85) 



and 



where Spx is a stochastic function of time that increments Px by ±Ap^ (equal chance for + 



or -) during each time step At, where Apx = ^j2D^^Kt and Dp^ is a constant. We work in 
the shock wave frame, with = —d^/dx, where 




We set Ey — —^sh^shBo/c, where Vgh, Ishi and Bq are the shock speed, the Lorentz factor 
for the shock speed, and the upstream magnetic field in the laboratory frame, respectively. 
Although Bz is in general not constant, we assume a constant value, B^ = 'jshBo- We 
set = 0.25c, '-fsh = 1.0328, and ^o/Bod = 4.5. The minimum of the approximate 
effective potential $ + xB^ occurs at x = 1.63d. We start the particles at a; = 2.25d, 
with Px = —O.lrriiC, Py = —lOrriiC, and Pz = 0. By varying Dp^, we are able to consider 
different values of a in equation (59). 

To improve the statistics, we implement particle splitting as follows. We start each 
simulation with 1000 particles. At each time step, we update p and x for each particle. 
When a particle escapes from the effective potential well the calculation for that particle is 
stopped and 7 at that time is recorded. When half of the original particles have escaped, 
each of the remaining particles is "split" — that is, a copy of each particle is created with the 
same position and momentum, and both the original and the copy are subsequently tracked. 
When we calculate the energy spectrum, the weightings of split particles are halved. The 
splitting process is repeated each time that the particle number drops to 500 until the end 
of the simulation, when every particle's 7 exceeds a threshold value, jm = 500. 

In the weak scattering case, a <C 1, a difficulty for calculating the spectrum is that 
almost all of the particles reach the threshold 7^ before escaping. We therefore use another 
method to obtain the spectrum when a <^ 1. We assume that the spectrum is given by a 
power law, dN/dj cx 7"'' with b > 1. The number of particles with 7 between 71 and 72 is 
given by 




(88) 



where C is a constant. The number of particles with 7 > 72 is 
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Prom Eqs. (88) and (89), we obtain the spectral index b as 



6 = 1 



ln(l + iVi^2/iV2^oo) 

In (71/72) 



(90) 



We take 71 = 250 and 72 = 500 and track 5000 particles for each value of a. 



Our numerical results for a < 1 are plotted in figure 5 along with our analytic result 
[equation (80)] for the same shock parameters. The numerical results and analytic result 
converge for a <^ 1. Our numerical results for a > 1 are plotted in figure 6, and are consistent 
with equation (84) when a ^ 1. 
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Fig. 5. — Power-law index b of the energy spectrum of accelerated particles as a function of 
the dimensionless momentum-diffusion parameter a. The crosses are the results of Monte 
Carlo calculations. The solid line is the analytic result in equation (80). 



4. Summary 



In this paper we model surfing acceleration of relativistic particles at perpendicular 
shocks as a process in which particles undergo diffusion in (x being the direction of the 
shock normal) while propagating in steady-state electric and magnetic field profiles. For 
shocks with > over some interval in x, where z is the direction of the magnetic 
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Fig. 6. — Log-log plot of the power-law index h of the energy spectrum of accelerated 
particles as a function of the dimensionless momentum-diffusion parameter a. The crosses 
are the results of Monte Carlo calculations. The dashed line gives the approximate result in 
equation (84). 

field, and for relativistic ions propagating primarily in the — y direction, surfing produces 
a power-law energy spectrum of accelerated particles, provided the momentum diffusion 
coefficient D^^ is independent of energy. We have calculated analytically the power-law 
index for the case of weak momentum diffusion, given in equation (80), and have carried 
out Monte Carlo test-particle calculations to determine the power-law index for the weak, 
moderate, and strong momentum diffusion cases (section 3). 
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